Distributional Ergodicity in Stored-Energy-Driven Levy Flight 
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We introduce a new class of random walk, the stored-energy-driven Levy flight (SEDLF), whose 
jump length is determined by a stored energy during a trapped state. The SEDLF is a novel 
type of anomalous diffusion, including the usual continuous-time random walk as a spacial case. 
We show that the ensemble-averaged mean square displacements exhibit subdiffusion as well as 
superdiffusion, depending on a parameter. We find that the diffusion coefficients of time-averaged 
mean square displacements are intrinsically random, a manifestation of distributional ergodicity. The 
diffusion coefficient shows aging in subdiffusive regime, whereas it increases with the measurement 
time in superdiffusive regime. 



Single particle tracking experiments in biological sys- 
tems often show that diffusion is not normal but rather 
anomalous [l|-|6|; that is, the mean square displacement 
(MSD) does not grow linearly with time but follows a 
power-law scaling 



(x?)oci^ (/3^1). 



(1) 



Because anomalous diffusions including subdiffusion 
(/3 < 1) as well as superdiffusion (/3 > 1) are ubiqui- 
tously observed in many biological experiments, anoma- 
lous diffusion is believed to play significant roles in cell 
biology such as gene regulation [7[ and active transports 
[l|, |6|. However, the underlying physical mechanisms re- 
main controversial. 

To understand the underlying mechanisms of these 
anomalous diffusions, phenomenological models such as 
continuous-time random walk (CTRW), Levy walk and 
flight, and other stochastic models of anomalous diffu- 
sion have been intensively studied [5|, IsMlSf . Among these 
models, CTRW shows a prominent feature called distri- 
butional ergodicity 10l4l2l. Il4|: that is, the time average 
of an observable converges to a random variable, i.e., con- 
vergence in distribution, but it does not coincide with the 
ensemble average as in the ordinary sense of ergodicity. 
It is considered that this distributional behavior of time- 
averaged observables in CTRW is related to large fluctu- 
ations of transport coefficients in single particle tracking 
experiments [2t5|. It is known that such distributional 
behavior is universal in infinite ergodic theory 



laui 



where ergodicity is satisfied with respect to an infinite 
(non-normalizable) invariant measure. This is a com- 
pletely different feature from other stochastic models of 
subdiffusion. 

Most of the studies on CTRWs so far are devoted to the 
case in which trapping time and jump length are mutu- 
ally independent. Effects of a coupling between them, 
however, become physically important for nonthermal 
systems such as cells [l|, l6J. In such nonthermal systems, 
a particle in a trapped state would not be simply frozen, 
but rather it would be storing a sort of energy for the 



next jump. Thus, a random walk driven by an stored en- 
ergy during a trapped state is essential in such nonther- 
mal systems, and it will also be important in complex 
systems such as finance 17[ and earthquakes 18 1. 

Here, we propose a class of CTRWs, which we refer 
to as the stored-energy-driven Levy flight (SEDLF), in 
which jump length and trapping time are mutually de- 
pendent. We show that the SEDLF exhibits novel types 
of anomalous diffusions and distributional ergodicity. In 
particular, it shows a whole spectrum of diffusion: sub- 
, normal-, and super-diffusion, depending on a param- 
eter 7, which characterizes the coupling strength be- 
tween jump length and trapping time. In addition, time- 
averaged observables such as the time-averaged MSDs 
(TAMSDs) are intrinsically random even when the mea- 
surement time goes to infinity. 

The SEDLF is based on CTRW with a non-separable 
joint probability of trapping time and jump length. In 
general, CTRW is defined through the joint probabil- 
ity density function (PDF) ■i/j{x,t), where ■i/j{x,t)dxdt is 
the probability that a random walker jumps with length 
[x, x-\-dx) just after it is trapped for period [t, t-\-dt) since 
its previous jump [l9|, |20[. In particular, the separable 
case il){x,t) = w{t)l{x), in which the jump length and 
the trapping time are mutually independent, has been 
extensively studied [8|, llCll4l2l |2l|. Here, we consider a 
non-separable case defined by 



ip{x,t) — wit) 



5(a;-tT) + (5(x-t-iT) 



(2) 



where ■w{t) is the PDF of trapping times and 7 e [0, 1] is 
a coupling strength. Note that a random walker under- 
goes a long trapped state before it performs a long jump. 
In addition, we assume that the PDF of trapping times 
follows a power law: 



w{t) 



Co 



(3) 



as i — > cxD. Here, a G (0, 1) is the stable index, a constant 
Co is defined by cq = c/|r(— q;)| with a scale factor c. 



For 7 = 0, the SEDLF is just a separable CTRW with 
jumps only to the nearest neighbor sites. On the other 
hand, for 7 > 0, the PDF of jump follows a power law: 



lix) 



'X '' 



i'{x,t)dt = 
27 



(n^) 



Co 



1 



27 |a;|i+"/T 



(4) 

Thus, the mean jump length diverges for 7 > a. Note 
that the Levy flight also has a power law distribution 
of jump length, which causes a divergence in the MSD. 
By contrast, the MSD of the SEDLF is finite with the 
aid of the coupling between jump lengths and trapping 
times as shown below. This property makes the SEDLF 
a physically more coherent model than Levy flight. 

Generalized master equations for CTRWs obtained in 
[19| can be utilized for our model. In general, the spaeial 
distribution P{x,t) of CTRWs with initial distribution 
Poix) at time zero satisfies the following equations: 

P{x,t)= f dt''i>{t-t')Q{x,t') + 'i>o{t)Po{x), (5) 



Q{x,t) 



dx' dt'ip{x',t')Q{x-x',t-t') 



dx'ipo{x\ t)Po{x — a;'). 



(6) 



where Q{x,t)dtdx is the probability of a random walker 
reaching an interval [x, x + dx) just in a period [t, t + di), 
^(i) is the probability of being trapped for longer than 
time t: 

/OO f^t /'t 

dx' / iix',t')dt' = 1- w{t')dt', (7) 
-00 Jo Jo 

^o{x,t) is the joint PDF for the first jump, and ^o(i) is 
the probability that the first jump does not occur until 
time t. Fourier-Laplace transform with respect to space 
and time (x — > fc and t — > s, respectively), deflned by 



gives 



P{k,s) 



P{k,s) 



ikx ^ — st 



dx / dtP{x,t)e"'''e 
, Jo 



Po(fc) l-wois)+(pQ{k,s) 



^{k,s) 



(8) 



(9) 



where (^0(^1 s) = [1 — w{s)]ipQ(k, s) — [1 — wo(s)]V'(fc, s). 
In the case of the SEDLF, we obtain i^^k, s) from 
Eq. (HI as follows 



ip{k,s) 



■ cos {kt'^) w{t)dt. 



(10) 



Note that -0(0,5) = w{s), and the asymptotic behavior 
of the Laplace transform of w{t) [Eq. ([3|)] is given by 



l~w{s)~cs°' (s^O). 



(11) 



We assume that the initial distribution Po{x) is the delta 
function, ^0(2;) = S{x), and w{t) = wo{t) (ordinary re- 
newal process [22|). As a result, we have the following 
generalized master equation in the Fourier and Laplace 
space: 



P(fc,s) 



1 l~w{s) 
s 1 ~7p{k,sy 



(12) 



where "0(^7 s) and w{s) are given by Eqs. (jlOp and (|11|) . 
Here, we derive the asymptotic behavior of the mo- 
ments of position xt for i — >■ cxd using the Fourier-Laplace 
transform P{k,s). The Laplace transform of (xt), de- 
noted by {xs), is given by (x^) ~ — i 

which means there is no drift, (xt) - 

the Laplace transform of the second moment, i.e., the 

ensemble- averaged MSD (EAMSD), is given by 



dP(k,: 



dk 

0. 



= 0, 

k=0 

Similarly, 



{-D 



d^P{k,s) 



dB 



fe=0 



Using the asymptotic behavior at s 
r(27-a) 1 



{^D 



\T{~a)\ 51+27 
|r(-a)|s"+i ^°^' 

„„a-|-l • 



The inverse Laplace transform for t 



1 0"(O,g) 
si — w(s) 

— !• 0. we have 

(27 > a) 

(27 = a) 

(27 < a) 
00 reads 



(13) 



(14) 



(X?) 



r(27 - a) 



|r(-a)|r(l + 27) 

1 

|r(-a)|r(l + a) 

\ / j-OL 

tcr(l + a) ' 



-e\ 



(27 > a) 
Plogi, (27 = a) (15) 
(27 < a) 



where we used {t^'^) = {P) = J^^x'^l{x)dx. Figure [T] 
shows the EAMSD for a = 0.5. Theory ^5^ is in excel- 
lent agreement with numerical simulations. 

It is also possible to derive higher order moments in 
the following way. By Eq. ((T2|) . we have the relation, 
P{k, s) = P{k, s)ip{k, s) + C, where C does not depend 
on k. Differentiating both sides n times with respect to 
fc, we have 



P(")(fc,s) 



ra-l 



^—-J2nClP^'\k,s)i;^--'Hk,s) 



(16) 




FIG. 1. Ensemble-averaged mean square displacements (q = 
0.5). Symbols are the results of numerical simulations for dif- 
ferent 7 with theoretical lines. There are no fitting parameters 
in the theoretical lines. We set the PDF of the trapping time 



as w{t) 



at 



(t > 1) in all the numerical simulations. 



Thus, the jump length PDF is given by l{x) = a/2^\x\^+"^'< 
from Eq. Q, and (l^) = a/ {a - 27) for 27 < a. 



From Eq. pUj) . we have -0(^"+^)(O, s) = 0. Accordingly, 
we obtain P^'^"+'^\0, s) = and 



^ n—1 

(17) 
by induction. Thus, (ccj""^^) = and the leading order 
for the Laplace transform of (xj") is given by 



i^T) 



M„(a,7) 

gl+2ri7 ' 

2"s[l -w{s)]" 



(27 > a) 



(27 < a) 



(18) 



where M„(a,7) is given by a recursion relation 
Mnian) = ELi2n^2iM„_,(a,7)r(2/7-a)/|r(-a)|. 
The above equations can be confirmed by Eq. (J17p and 
mathematical induction. Therefore, the asymptotic be- 
havior for s — >■ is given by 



(^r) 



Af„(a,7) 

^l+2n7 ' 

(2n)!{log(^)}" 

2"|r(-a)|"s""+i' 

I (2c)"s""+i ' 



(27 > a) 

(27 = a) (19) 

(27 < a) 



and the inverse Laplace transform for i — > 00 reads 



(-?") 



M.(«,7) 



r(i 



,/2n7 



2^7) 

(2n)! 



2"|r(-a)|"r(l- 
[ (2c)"r(l-Hna)' 



no) 



<""{logi}". 



(27 > a) 
(27 = «) 
(27 < a) 



(20) 

It follows that the distribution of a scaled position 
Xt/ yjxf) converges to a time-independent non-trivial 
distribution. In other words, Xt/ ^ {xf) converges in dis- 

^a,7 

' °° (zfc)2»M„(a,7)r(l + 27)^ 



tribution to Xa^^ as t — ;■ 00, where 



ikXc, 



E 



'^> = 



(2n)!Mi(a,7)"r(l + 2?i7)' 



^ (zfc)^"r(l -f g)-^ 

^ 2"r(i 



no) 



(27 > a) 



(27 < a) 



We note that the distribution of the random variable 
Xa.-y for 27 < a is called a symmetric Mittag-Leffler 
distribution of order a/2 [23|, |2J]. Figure [5] shows the 



PDFs of Xa^^ for 7 = 0.1,0.6, and 0.8 {a ^ 0.5). The 
PDF depends crucially on 7 when 27 > a. 

Here, we investigate ergodic properties of time- 
averaged MSD (TAMSD), defined by 



<52(A;f) 



1 



t-A 



t-A 



[x{t' + A)-x{t')]^dt'. (22) 



It has been known that TAMSD can be represented using 
the total number of jumps [1^ l24|, i.e., Nt, and hk = 



^k-l 



^'fe + 2 J2m=l lklmS{A — tk+ tm)- 

1 ^* 

52(A;t)^-^/ife (i^^), (23) 

fc=0 

where Ik is the fcth jump, tk is the time fcth jump oc- 
curs, and 9{t) is a step function, defined by 6{t) = for 
t < and t otherwise. One can show that X]fc=o(^'t ~ 
All)/ J2klo ll -^ aa t ^ 00. It follows 



(52(A; t) ~ DtA (A < t and t ^ 00) 



(24) 



where A = Eflo^fc/*- ^s shown in Fig. El TAMSDs 
increase linearly with time A (normal diffusion), while 
the diffusion coefficients show large fluctuations. 

Now, we derive the PDF P2{z,t) of Zf = 'Ekio^l- 
We note that If and Nt are mutually correlated because 
both of them depend on the fcth trapping time, and thus 
we cannot apply the method used in previous studies 



12l |24| . Instead, we use the fact that Zt obeys a directed 

(25) 



SEDLF with the joint probability 
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FIG. 2. Probability density functions of a scaled position Xt = xtj \J {xl) (q = 0.5). PDF P(Xt) converges to a non-trivial 
PDF as t — > cxD. (a) P{Xt) converges to a symmetric Mittag-LefSer distribution with a = 0.5 (7 = 0.1). (b) and (c) The PDFs 
of Xt converge to different distributions. The PDF wify used in the numerical simulation is the same as that in Fig. [l] 




FIG. 3. Time-averaged mean square displacements {a — 0.5 
and t = 10 ). TAMSDs for eight different realizations are 
drawn in (a), (b), and (c) for 7 = 0.1,0.6 and 7 = 0.8, re- 
spectively. Linear scalings are shown by the solid lines for 
reference. The PDF w{t) used in the numerical simulation is 
the same as that in Fig. [T] 



Therefore, in the same way as Eq. (|12|) . we obtain 



A(fc,.) = i-'^^(^) 



Sl-V2(fc,s) 



(26) 



with ip2{k,s) = /q°° e~'**e''''* '' w{t)dt. Thus, the calcula- 
tions of the moments Zt are almost parallel with the case 
of Xt- For example, we obtain the mean diffusion coeffi- 
cient, {Dt) — {Zt)/t, as follows: 



r(27-a) ^21- 



|r(-Q)|r(n-27) 



(27 > a) 



(A) 



\T{-a)\T(l+a) 



cr(i-i-Q) 



r-ilogf, (27 = a) (27) 

(27 < a) 



for i ^' 00. For 27 > 1, the diffusion coefficient en- 
hances, otherwise it shows aging. This enhancement of 
diffusion coefffcients is completely different from separa- 
ble CTRWs [13, El Q H and correlated CTRWs [H], 
both of which show only aging. 



Furthermore, the second moment of Dt is given by 



(A^) 



r(47-^a)|r(-a)|+2r(27-a)^ ,4T-2 
r(47+l)|r(-a)|2 *- 



r(2Q+i) 

yr 

c^r(2Q+i; 



2a 



(27 > 


a) 


(27 = 


a) 


(27 < 


a) 




(28) 



In a way similar to the calculation for xt, we can ob- 
tain all the higher moments of Dt- Therefore, the 
scaled diffusion coefficient Dt/yiDf) converges to a 
time-independent non-trivial distribution as t —>■ 00, in- 
dicating that the scaled diffusion coefficient converges 
to a random variable (i.e., distributional ergodicity). In 
fact, the relative standard deviation (RSD) of Dt, cteb = 
\/(D?) — (Df)'^ I (Df) , which is an ergodicity breaking pa- 
rameter [ll|, |l2|, |lj, |26|, |27| , remains constant as i — > 00: 



C^EB 



r r(47-a)|r(-a)| ^ 
[ r(27-a)- +^ 



V2*(a) - 1 



$(27) - 1 (27 > a) 

(27 < a), 
(29) 



where $(a;) = r(.T + if /T{2x + 1). As shown in Fig. IH 
the RSDs of Dt depend on 7, and they are different from 
that in CTRW when 27 > a. 

In conclusion, we have provided a novel random walk 
that generates subdiffusion as well as superdiffusion. 
This model (SEDLF) removes unphysical situations in 
Levy flight such that the EAMSD always diverges. In 
the SEDLF, we have shown that time-averaged diffusion 
coefficients converge in distribution (distributional ergod- 
icity) . Distributions of the diffusion coefficients depends 
not only on the exponent a of the trapping-time distri- 
bution but also on the coupling exponent 7 for 27 > a, 
and thus are different from those in separable CTRWs 
111 ] as well as random walks with static disorder 114 1. 



Especially, in superdiffusive regime (7 > 0.5), the mean 




Y 

FIG. 4. Relative standard deviation of Dt as a function of 
J {a — 0.25,0.5, and 0.75). Symbols are results of numeri- 
cal simulations. We calculate Dt by 5^{A;t)/A in numerical 
simulations with A = 10 . Solid lines are theoretical curves 
(|29[) . The PDF it)(t) used in the numerical simulation is the 
same as that in Fig. [1] 



diffusion coefficient enhances according to the increase of 
the measurement time. 
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